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ABSTRACT 

We investigate the variability timescales in the jet of M87 with two goals. The first is to 
use the rise times and decay times in the radio, ultraviolet and X-ray lightcurves of HST-1 
to constrain the source size and the energy loss mechanisms affecting the relativistic electron 
distributions. HST-1 is the first jet knot clearly resolved from the nuclear emission by Chandra 
and is the site of the huge flare of 2005. We find clear evidence for a frequency-dependent decrease 
in the synchrotron flux being consistent with energy losses. Assuming that this behavior is 
predominantly caused by synchrotron cooling, we estimate a value of 0.6 mG for the average 
magnetic field strength of the HST-1 emission region, a value consistent with previous estimates 
of the equipartition field. In the process of analyzing the first derivative of the X-ray light curve of 
HST-1, we discovered a quasi-periodic oscillation which was most obvious in 2003 and 2004 prior 
to the major flare in 2005. The four cycles observed have a period of order 6 months. The second 
goal is to search for evidence of differences between the X-ray variability timescales of HST-1 and 
the unresolved nuclear region (diameter < 0.6"). These features, separated by more than 60 
pc, are the two chief contenders for the origin of the TeV variable emissions observed by HESS 
in 2005 and by MAGIC and VERITAS in 2008. The X-ray variability of the nucleus appears to 
be at least twice as rapid as that of the HST-1 knot. However, the shortest nuclear variability 
timescale we can measure from the Chandra data (< 20 days) is still significantly longer than the 
shortest TeV variability of M 87 reported by the HESS and MAGIC telescopes (1 - 2 days). 

Subject headings: Galaxies: active — galaxies: jets — galaxies: individual (M87) — X-rays: general — 
radio continuum: galaxies — radiation mechanisms: nonthermal 
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1. Introduction 

We have been monitoring the jet of M87 in the 
X-rays with Chandra and the ultraviolet (UV, at 
A=220 nm) with HST since 2002 Jan, and in the 
radio since 2003 {VLA, primarily at 15 GHz) and 
2005 ( VLBA, primarily at 1.7 GHz). Previous pa- 
pers from this project include Paper I reporting 
our first results (Harris et al. 2003), Paper II 
which focused on the HST data (Perlman et al. 
2003), Paper HI which was mainly on the X-ray 
lightcurve of HST-1 which delineated the massive 
flare in 2005 (Harris et al. 2006) and Paper IV, the 
VLBA results showing superluminal proper mo- 
tions in HST-1 (Cheung, Harris & Stawarz 2007). 

In this paper (V of the series), we present an 
analysis of the lightcurves for the nucleus, the jet 
knot HST-1 which is 0.86" from the core and which 
is the site of the massive X-ray, UV, and radio fiare 
described in Paper HI, knot D, and knot A. Knot 
A was mainly used as a control source since we 
do not expect short timescale variability because 
it is well resolved in all bands. The lightcurve of 
Knot D illuminates the effects of HST-1 on adja- 
cent regions since it appears to have very little if 
any intrinsic variability on the timescales of inter- 
est here. In § [2l [3l and [4] we describe the data 
and the analyses methods. In §[5] we use the rising 
segments of the lightcurves to derive upper limits 
on the size of the emitting regions, and in § [5] we 
examine the decay timescales of HST-1 in X-rays, 
UV, and radio bands in order to isolate signatures 
of E'^ losses. We describe newly discovered oscilla- 
tions in the brightening and fading of HST-1 in §[71 
Finally, in § [8l we discuss the evidence for short 
timescales in the X-ray variability of the nucleus 
and HST-1 which is relevant to the question of the 
location of the variable T eV emission reported b y 



the HES S collaboration ( Aharonian et al. 2006h . 



MAGIC (I Albert et all EoOSh. and VERITAS (e.g 



Ergin 20081 ). Some p reliminary results fr om this 
work were reported in [Harris et al. I (|2008h . 

We take t he distance to M87 to be 16 Mpc 
(jTonrvl Il99lh so that 1" corresponds to 77 pc. 
Throughout this paper we assume that the ra- 
dio to X-ray emission from all parts of the M87 
jet comes from synchrotron emission, as argued in 
our previous papers (I & III). In particular, we 



^Astronomical Observatory, Jagiellonian University, ul. 
Orla 171, 30-244 Krakow, Poland. 



assume the X-ray nuclear emission is dominated 
by synchrotron emission from components of the 
inner (unresolved) jet rather than by thermal pro- 
cesses associated with the accretion disk around 
the central black hole. 

2. The Chandra X-ray Data 

Since 2002, we have observed M87 with Chan- 
dra 6-7 times each observing season with 5 ks ex- 
posures typically separated by six weeks. Addi- 
tionally in 2005, near the maximum of the X-ray 
lightcurve of HST-1, we scheduled weekly obser- 
vations to constrain shorter timescale variability 
(epochs Ya to Yg; see Tabled]). After the report 
of variable TeV e mission in the 200 5 HES S obser- 
vations of M87 bvlXh aronian et al.l (l2006l ). we ob- 
tained Director's Discretionary Time (DDT) ob- 
servations to sample the X-ray lightcurves on ~2- 
3 day timescales during two 'dark time' fortnights 
in Feb. and Mar. of 2007 when TeV observations 
were scheduled (epochs Ys to Zb). A total of 61 
observations have been obtained from these pro- 
grams thus far. 

For details of our reduction procedures, see Pa- 
pers I & III. Briefly, we use a l/8th segment of 
the back illuminated S3 chip of the ACIS detec- 
tor aboard Chandra. This permits us to have 
a frame time of 0.4s with 90% efficiency. Al- 
thoug h this setup wa s essen tially free of pileup 
when IWilson fc Ya^ (|2002l ) tested various op- 
tions during 2000 July, with the advent of the 
ever increasing brightness of HST-1, pileup ( David 
20011 ) became a major problem so we switched to a 
detector based measure of intensity: keV/s. This 
approach uses the event 1 file with no grade filter- 
ing (so as to recover all events affected by 'grade 
migration') and we integrated the energy from 0.2 
to 17 keV so as to recover all the energy of the piled 
events. Other uncertainties for piled events comes 
from the on-board filtering, the 'eat-thy-neighbor' 
effect, and second order effects such as release of 
trapped charge (see § [X| . 

2.1. Photometry 

Although we used small circular apertures for 
fluxmap photometry in Paper I, the basic analy- 
sis for this paper adopts the rectangular regions 
used in Paper III so as to encompass more of the 
point spread function (PSF). The 4 regions of in- 
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terest (the core, HST-1, knots D and A) are shown 
in fig. m we did not use background subtraction 
because the photometric apertures were small. 

All events within each rectangle are weighted 
by their energy and the sum of these energies, 
when divided by the exposure times, gives the fi- 
nal keV/s value used in the lightcurve. Uncertain- 
ties are strictly statistical, based on the number 
of counts measured: ^fN /N and typically range 
from 1% to 5%. 

To analyze and compare timescales, we measure 
the slope between adjacent measurements (and 
also between every other observation) by calcu- 
lating the ratio, {I2 — Ii)/At, where Ii and I2 
are the intensities at the times ti and t2 and 
At[yrs] ~ t2 — ti. To convert this to a fractional 
change we divide by min(/i,/2), so the definition 
of fpy is: 



fpy 



D 

TKt' 



(1) 



where D = {I2 — Ii)- When the intensity is in- 
creasing {I I < I2), the fractional change per year 
has i=l and is denoted (when required), as fpy~^. 
When the intensity is dropping (/i > I2), i=2 and 
we specify by using fpy~. 

(Hereafter, the superscripts are suppressed where 
the signs are implicit). According to this defini- 
tion, the doubling time for a given value of fpy is 
simply tdoubie [yr] = ^/ fpy, thus when fpy = ±1, 
there was a rate of change which would produce a 
factor of two increase or decrease in / in one year. 

The uncertainties in each value of / are prop- 
agated to the first derivative by calculating the 
square root of the sum of the squares of the errors 
on intensity. Denoting cTj as the uncertainty of Ii, 
we express the error of fpy as: 



<^ifpy) 



—— X — X 
A< Ii 



Cj-Y + CfY (2) 



Here i=l, ']—2 for the error on fpy^, and i=2, 
for fpy~. 

The fpy values are plotted in fig. [5] for the nu- 
cleus, HST-1, knots D and A. The latter two fea- 
tures serve as controls, and display the expected 
behavior for a steady source: at long time inter- 
vals, fpy values are consistent with zero and have 
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Fig. 1. — A Chandra X-ray image constructed 
from the sum of our data from 2002 to 2004. This 
is a total power image for which each event has 
been multiplied by its energy. In the top panel, 
a Gaussian smoothing function of FWHM= 0.25" 
has been applied. Contours start at 0.01 eV s^^ 
per 0.049" pixel and increase by factors of two. 
The second panel shows the regions used for pho- 
tometry. The four rectangular regions in red (from 
left to right) are the nucleus, the flaring knot HST- 
1 (0.86" from the nucleus), knot D, and knot A. 
The long thin cyan rectangles and the dotted ma- 
genta rectangles are used for 'readout streak pho- 
tometry' ('on' and 'background', respectively) - 
see § 12.21 The image is a 5 ks exposure from 2005 
May when HST-1 was close to its peak intensity. 
Pixel randomization has been removed and in this 
figure panel, the events are binned into 0.123" pix- 
els. The axes are J2000 coordinates. 
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small uncertainties. The errors increase as At de- 
creases. Since all errors are la (^/N/N), some 
values for knots D and A appear to be different 
from zero, as expected. 

In the next section we discuss using the read- 
out streak for estimating source intensity and in 
the Appendix {§E| we describe some of the other 
problems engendered by pileup. 

2.2. Photometry of the Read-out Streak 

The only check we have been able to de- 
vise is 'read-out st reak photometry' (see also 



Marshall erallbOOSf ). For this procedure, we iso- 
late the segments of the readout streak which 
are not close to the jet. With long thin rectan- 
gles (nx2 pixels; fig. H]), and adjacent rectangles 
to measure the background, we estimate the effec- 
tive exposure time for the net counts as (^frames) 
X n X 41/is, where n=length of rectangles in pix- 
els. In a typical 5 ks observation, we are thus able 
to get the equivalent of 30 to 40 seconds worth of 
continuous clocking mode data which are free of 
pileup effects but suffer from poor s/n because of 
the necessarily large background area. While the 
general behavior of the streak photometry is con- 
sistent with the expectation that it would become 
increasingly larger than the standard photome- 
try as the intensity of HST-1 increased (more 
pileup and more on-board rejection), there are 
a few unexpected departures from this behavior. 
In particular, there are some high values where 
both methods give the same intensity, and there 
is one segment of low intensity where the streak 
photometry is significantly less than the standard 
photometry. For these reasons we rely on our 
'standard' photometry and relegate the results of 
the readout streak photometry to the status of 
'caveats' and 'alternate possibilities'. 

3. The UV Data 

The UV data used in this paper were obtained 
from a series of .ffST proposals (Biretta, PI) which 
were synchronized with the Chandra observations 
although for various reasons not every Chandra 
observations has a corresponding HST observa- 
tion. The data were reduced with the usual pro- 
cedures and will be described in a separate paper 
(Biretta, in preparation). 
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Fig. 2. — The first derivative of the lightcurves for 
the nucleus, HST-1, knots D and A (top to bot- 
tom). The smaller squares are values from inten- 
sity pairs which skip the intervening observation 
and are thus not independent of the larger circles 
(adjacent observations). In the upper panel the x 
axis is the time between observations; in the lower 
panel it is the date. The y axis is fractional change 
per year. Note the change in y scale for different 
plots in the upper panel, but in the lower panel 
we have set the max/min values to ± 15. In the 
lower panel, intervals with closely spaced observa- 
tions have been averaged to avoid large error bars 
arising from small At. 



4. The Radio Data 

The VLj^ 15 GHz observations were obtained 
as part of our multi-frequency program coordi- 
nated with the Chandra and HST monitoring and 
began mid-2003. In each VLA cycle, we observed 
M87 in three 8 hr runs (program codes AH822, 
AH862, AH885, AC843): two in A array followed 
by one in B array. The beam sizes were typically 
around 0.15" and 0.4", respectively. The longer 
gaps occur during C and D array configurations 
when the angular resolution is not sufficient to sep- 
arate the nucleus and HST-1. 

The observations utilized two adjacent 50 MHz 
wide intermediate frequencies centered at 14.94 
GHz. A total of 1-1.5 hrs of on-source time was 
split between 9-10 scans over each 8 hr run to 
obtain good (m, i;) coverage. The data were cali- 
brated in AIPS (Brid le fc Greisenl ll994^ with the 
flux density scale set using scans of 3C286 and the 
initial phase corrections determined with a nearby 
calibrator. Subsequent phase and amplitude 
self-calibration was performed using the Caltech 
DIFM AP package (jShepherd. Pearson. &: Taylor 
19941 ). 

An additional archival B-array dataset (BK073) 
was analyzed to give us the data point in January 
2000, before our monitoring began. The observa- 
tion used an identical setup to ours but with only 
four 3.5 min scans obtained. This was sufficient 
to detect a faint (2.8 mJy) feature at the position 
of HST-1, thus providing a baseline measurement 
to our subsequent ones. 

5. Estimating the Size of Emitting Vol- 
umes from Rise Times 

As described in earlier papers of this series, the 
rise of the light curve can give estimates of the 
source size so long as the beaming parameters - 
the bulk Lorentz factor of the jet (F), the angle 
between the jet and our line of sight (6), and the 
Doppler beaming factor (S) - are not changing. 

For the standard analysis of the X-ray light 
curve of HST-1, the largest observed slope oc- 
curred early in 2004 and had a value of fpy = 8 
(fig. H]); the intensity doubled during our 6-week 



sampling interval. Thus we are able to set an up- 
per limit to the characteristic size of the emitting 
region in the jet frame of diameter < 0.12(5 light 
years (45 light days). Significantly larger values 
of fpy occur for the HST-1 light curve generated 
from the readout streak photometry (§ 12. 2[) . but 
these have large uncertainties so do not yield use- 
ful limits on the source size. 

The smallest directly measured value for the 
size of HST-1 comes from the VLB A beamsize of 
3 mas which corresponds to 0.7 light years (Paper 
IV) . If the X-ray and radio emitting volumes were 
to be one and the same (i.e., the upstream end of 
HST-1), and if the radio size is actually similar to 
our VLB A resolution, then 5 would be of order 5, 
a value similar to that estimated in our previous 
papers. 

Because of the second order effects which were 
prevalent when HST-1 was strong, we have re- 
stricted our analyses of the nuclear fpy data to 
(a) 2000-2003 plus 2008 when HST-1 was at a low 
intensity (less than 2 keV/s), and (b) 2006-2007 
when the HST-1 intensity was at an intermedi- 
ate value (2 to 4 keV/s). The data between 2 
and 4 keV/s are used because this time interval 
includes the closely spaced observations of 2007. 
We believe the second order effects for this inten- 
sity regime are minimal because although the light 
curve of the nucleus shows a very obvious featue 
corresponding to the peak of HST-1 in 2005, there 
is no apparent correlation between the intensity of 
the nucleus and that of HST-1 for the 2006-2007 
data. These data are shown in fig. [3l 

During low intensity intervals, the maximum 
value of fpy^ was about 12 (a light travel time 
of 30(5 light days), although there were no closely 
spaced observations during these periods so that 
larger fpy^ values may have been missed because 
of inadequate sampling. When HST-1 was at in- 
termediate levels (~2-4 keV/s), we had the closely 
spaced observations during 2007 Feb and Mar from 
the DDT project and the maximum fpy^ for the 
core was 28.5±5.6, but we take a value of 19 (light 
travel time: 196 light days) as a characteristic 
value since the second largest value was 18.8±1.1. 



^The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under coopera- 
tive agreement by Associated Universities, Inc. 
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Fig. 3. — fpy values for the nucleus. The 
larger open circles come from adjacent observa- 
tions whilst the smaller squares are data pairs 
formed by skipping the adjacent observation. First 
panel: the period from 2000-2004 and 2008 when 
the HST-I intensity was low (<2 keV/s). Second 
panel: the period 2006-2007 when the intensity of 
HST-1 was between 2 and 4 keV/s. 
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Fig. 4. — The X-ray, UV, and radio lightcurves of 
HST-1. The intensity is plotted on a log scale 
to demonstrate the overall conformity between 
bands. Each curve has been normalized by setting 
the peak value to unity so as to permit visual com- 
parison of the decay. Peak values are 12.417 keV/s 
(X-ray); 0.596 mJy (UV); and 0.084 Jy (15GHz). 

6. Analysis of the Decay Phases of the 
Lightcurves of HST-1 

In this section, we analyze the decay of the light 
curves of HST-1 at X-ray, UV, and radio wave- 
lengths (fig. HI). We will not attempt to make a 
parallel investigation of the nuclear emission be- 
cause we have no information as to the size of the 
emitting volume or its geometry, the interpreta- 
tion of the UV data would be problematic, and al- 
though likely, it remains to be demonstrated that 
the nuclear X-ray emission is actually non-thermal 
emission from the inner jet rather than from some 
thermal process associated with the accretion disk 
or its environs. 

The decay of the light curve may be caused by 
several effects but has the potential to reveal which 
processes are dominant. If all bands have similar 
rates of decreasing intensity, the most likely cause 
is either a change in the beaming factor (as might 
arise from a change in 9) or a general expansion 
which reduces the energy of all electrons accord- 
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ing to their energy (the so called "E^ losses") as 
well as reducing the magnetic field strength. Note 
however that if the emitted spectrum is not a sim- 
ple power law, but for example steepens at high 
frequencies, then a simple expansion may produce 
a much stronger decay at high frequencies both be- 
cause the previously 'viewed' electrons for a fixed 
observing band are now replaced by the fewer (pre- 
viously higher energy) electrons and also because 
of the weaker B field, the fixed observing band now 
comes from an even higher energy segment of the 
electron distribution. 

Our preliminary analysis of Paper III indicated 
that the initial decrease of the major flare had a 
similar timescale as the preceding rise, and that 
this might be indicative of either a changing 6, or a 
compression and subsequent expansion. However, 
it has become clear that although the UV and X- 
ray light curves appeared to have a similar behav- 
ior initially, there are instrumental effects present 
which were not recognized, there are significant 
differences in the UV and X-ray decays and the 
radio intensity did not conform to the rapid decay 
seen at higher frequencies (fig. H]). 

We investigate two aspects of this problem. 
First, we examine the behavior of the UV and ra- 
dio lightcurves at times when large rates of decay 
are observed at X-rays, and second we compare the 
fpy values between bands without regard to when 
they occurred. Both of these approaches suffer 
from the sparser sampling in the radio. Although 
most of the UV data were obtained within a week 
of the Chandra observations, not every Chandra 
observation has a corresponding TJST observation. 

6.1. Comparison of particular time seg- 
ments 

This approach is based on the assumption that 
all 3 wavelength bands come from the same emit- 
ting volume. If however, the emitting volume 
is 'layered' (e.g., concentric spheres with longer 
wavelengths coming from larger volumes), then it 
would be possible to have different characteristic 
decay times for each layer (assuming different val- 
ues of the magnetic field strength), and the de- 
cays in the light curves would not have to happen 
at the same time since the cessation of injection 
of particles and fields would not necessarily be si- 
multaneous in all emitting volumes. 



The most significant (i.e., without excessively 
large uncertainties) decays in the X-ray lightcurve 
occurred in (a) 2002.0 {fpy = -3.32 ± 0.45), 
(b) 2005.5 {fpy = -3.51 ± 0.08), and (c) 2007.7 
(/py = -4.53 ±0.85). 

For event (a), we have no relevant radio values 
of fpy and the UV sampling was not sufficient: 
the UV intensity rose from 28 to 52 /iJy between 
2001.58 and 2002.16. The decay in the X-rays oc- 
curred between 2002.044 and 2002.119. 

Event (b) occurred after the peak of the giant 
flare. This is perhaps the best example of what 
might be expected from i?^ losses affecting the 
electron energy distribution (see fig.|l]). However, 
as mentioned above, we can't rule out the possi- 
bility of an expansion of the source if the X-rays 
are coming from a segment of the electron distri- 
bution that is falling more rapidly than the power 
law connecting the UV to the X-ray. The fpy 
value for the UV is -0.85±0.86. 

The third notable decline in the X-ray light 
curve came during the 3 months in 2007 whilst 
M87 was behind the Sun (Aug. - Nov.). Unfortu- 
nately, the UV monitoring had a much longer gap 
(almost 7 months) so the relevant data for a direct 
comparison do not exist. 

If the X-ray fpy~ value of -3.5 during the main 
decay phase ('event b', fig. [5]) was dominated by 
E'^ losses, we can make an order of magnitude es- 
timate for the magnetic field strength. Note that 
we do not consider here a detailed evolution of the 
synchrotron flux produced by the radiatively cool- 
ing electron ene rgy distribution as considered in, 
e.g., iKardashevI (_1962), but we can approximate 
the cooling time required to drop the intensity by a 
factor of two by assuming that this decay is caused 
by a factor of two fewer electrons with energies 
providing the bulk of the observed X-rays between 
0.2 and 6 keV. Further, we assume the exponent 
of the electron distribution, p=2aa;-l-l=3.4, which 
is our best estim ate for the spectru m of HST-1 in 
the X-ray band (|Harris et al.ll2006f ). Since dE/dt 
to first approximation shifts the power law dis- 
tribution, N(E) to lower energies, we require an 
energy shift of 2~^^^'^ or 0.8. We then ask what 
magnetic field strength, B is required to produce 
a 20% energy loss in the time, r, it takes for the 
intensity to fall by a factor of two. The observed 
time. To is l/3.5=0.28yr; in the jet frame t'=6to 
where S is the beaming factor. 
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Fig. 5.— Multiband data for HST-1 during 2005. 
The top panel shows the Ughtcurves: circles for 
X-ray; squares for UV with dashed line; and tri- 
angles for the 15 GHz data (dotted line). The 
lower panel shows the corresponding fpy values, 
the first derivative of the light curves. Coding is 
the same as the upper panel except no line is added 
to the 2 radio data, both of which are consistent 
with zero. Note the shoulder in the UV lightcurve 
following the peak. This can also be seen in fig.IH 
and leads to the much smaller values of dl/dt for 
the UV than for the X-ray. The radio lightcurve 
is essentially 'flat topped' and shows no change 
during the year. 

Wi th the stand ard e quations (3.28 and 3.32) 
from [Pacholczyls (|l97Gf ). setting r'dE/dt=0.2E 
(where E is the energy of the electron) , and chang- 
ing jet frame parameters to the observer frame, we 
find 



B « 0.5^1/3 ( 



-2/3 



(kev) 



-1/3 



iG . (3) 



Here To/yr = 1/ fpy^ is the observed time for the 
intensity to drop by a factor of two and e is the 
characteristic energy of the X-ray band. For our 
parameters, this reduces to B(5^/"^=l.l mG, and for 
(5=5, 0.6 mG, a value reasonably consistent with 



the 1 mG derived on the basis of equipartition 
conditions for HST-1 before the major flare (Paper 

I)- 

If losses were the controlling factor in the 
light curve decay, we would expect r(UV) to be 
longer by the square root of the ratio of the fre- 
quencies, «14. The actual observed t(UV) is 
1/0.85 = 1.18 yr although the uncertainties in- 
clude much longer times {fpy^O is within the la 
error). Of course the fact that both the UV and 
the radio intensities decline significantly at later 
times indicates that expansion of the source is 
probably a major contributor to the light curves 
behavior during some intervals. 

6.2. Comparison of extreme values of dl/dt 

In fig. [S] we plot the fpy values for HST-1 in 
different bands. Although many of the uncertain- 
ties are large, and the sampling in the radio is 
clearly insufficient, we find that the largest (ab- 
solute value) believable negative fpy's are -5, -2, 
and -1 for the X-ray, UV, and radio, respectively. 
We take this as evidence that energy loss by expan- 
sion alone, is not indicated. The obvious caveat to 
this conclusion is that there is a spectral break in 
the optical/UV in the sense that > ctro- If 
there is additionally a curving downwards of the 
spectrum in the v — 10^^ to 10^^ Hz band, we 
could explain the fpy data with just expansion. 

7. Discovery of Impulsive Brightening in 
HST-1 via dl/dt 

Close examination of figs. [2]and[7] (upper panel) 
shows that the X-ray light curve has a series of 
small peaks superposed on a gradually rising in- 
tensity between 2002 and 2004. Their reality is 
demonstrated by the plot of the first derivative 
(figs. [2] & [7]). Characteristic times of this oscil- 
lation (peak-to-peak or trough-to-trough) range 
from 0.50 years (most common) to a maximum 
value of 0.84 year. Although there are some 
slightly discernible features on the /py(UV) plots, 
the data are not sufficiently numerous or robust 
enough to search for lags between bands. The 
impulses are not evident at radio frequencies and 
their existence is debatable for the UV although if 
they were as large in the UV as in the X-ray, they 
should have been detected. 

Although the causes of these oscillations are not 
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Fig. 6. — The X-ray, UV, and radio scatter plots of 
fpy values for HST-1. The smaller square points 
are for every other observation and are thus not 
independent of the data for adjacent observations 
(larger circles). Top to bottom: X-ray, UV, and 
15 GHz. The absence of large values at large times 
in the top panel arises from the smaller intervals 
between observations: if we had calculated fpy for 
every possible pair of observations, we would have 
recovered the data describing the giant flare. 



known, we speculate that a quasi periodic varia- 
tion in the conversion of bulk kinetic jet power 
to the internal energy of the radiating plasma is 
more likely than a modulation of power flowing 
down the jet. We also disfavor a changing beam- 
ing factor caused by a changing angle to the line 
of sight (e.g. a "thrashing jet"). For jet modu- 
lation, a thrashing jet, and periodic compression 
and expansion, we would expect any oscillation to 
be evident at all frequencies equally. The absence 
of oscillations in the UV encourages us to look for 
an oscillating injection of particles which, for the 
highest energy electrons is made manifest by the 
short lifetimes, but for the electrons radiating at 
lower frequencies, gets smoothed out by the con- 
tinuously accumulating total number of radiating 
particles. 



Fig. 7. — Upper panel: lightcurves of HST-1: 
X-ray (solid), UV (dashed) and radio (dotted). 
Lower panel: fpy values for HST-1 prior to the gi- 
ant flare. Larger circles are the primary data from 
adjacent observations. Smaller circles are gener- 
ated from every other observation. The curves 
connect primary data. 

8. Comparing X-ray timescales of HST-1 
and the nucleus to evaluate the site of 
TeV flaring 

To date, there have been two reports of TeV 
short-timescale flaring f r om M 87: 2005 Apr 



(HESS, lAharonian et al 



(MAGIC, 



Albert et all 



20061) and 2008 Feb 
2008h : (VERITAS, e.g. 



lErgi ^ l2008i r While the angular resolution of the 
TeV systems does not permit a location to be 
determined, the expected relation between X-ray 
intensity and TeV intensity via an inverse Comp- 
ton model holds the potential of localizing the site 
of the TeV emission if an unambiguous feature in 
the TeV light curve can be associated with one in 
the X-rays. Moreover, since both instances of TeV 
flaring appeared to be characterized by timescales 
of only a few days, it is also possible to evaluate 
statistical differences in X-ray timescales, partic- 
ularly for the two leading contenders, the nucleus 
and HST-1. 

There are a few striking differences in fig.[2]be- 
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Fig. 8. — The first derivative of the nucleus from 
2006 during which time HST-1 was <4 keV/s. The 
top panel is a segment of the light curve of HST- 
1 and the middle panel shows the light curve for 
the nucleus with 5% of HST-1 subtracted. There 
are no egregious signatures of contamination. The 
bottom panel shows the first derivative of the nu- 
clear lightcurve with smaller squares for intensity 
pairs with an intervening observation. 



tween the nucleus and HST-1. For the nucleus, 
there are quite large values of fpy at short time 
scales, whereas HST-1 has an order of magnitude 
smaller amplitudes and these occur at somewhat 
longer timescales. We interpret the presence of 
large amplitudes at short sampling times together 
with smaller amplitudes at longer times to mean 
that we can characterize the nuclear variability as 
a sort of 'fiickering'. HST-1 of course provided 
us with a major flare with a timescale of a year 
or more (figs. I4|9p . The energy emitted by the 
nuclear fiickering is a small fraction of the energy 
emitted by the fiaring of HST-1, but the timescales 
are quite different. 

In Paper IV, we argued that the site of the 
TeV fiaring observed by HESS could be HST 



1, wh er eas others [e.g. iGeoreanopoulos et al 



(120051). [Ta7ecchio fc Ghisellinil (|2008h . and 



see 



also Levinsonl (|200Cll )] have suggested a location 



closer to the super massive black hole (SMBH). 
We pointed out that besides the coincidence of the 
peak of the HST-1 (X-ray, UV, radio) light curve 
occurring at the same time as the HESS event in 
2005, we knew that HST-1 was physically small 
(from the VLBA observations and from the X-ray 
variability), that the emitted power at TeV ener- 
gies was comparable to the X-ray power (linked, 
for example, by an SSC model), and that there 
were difficulties of getting TeV photons out from 
the immediate vicinity of the SMBH. While none 
of these considerations have changed, our cur- 
rent analyses on timescales could be thought of 
as circumstantial evidence that the nuclear X-ray 
emission comes from a smaller emitting volume 
than that of HST-1 and that this smaller emission 
region could be the same as that producing the 
TeV fiares, even though our upper limits on the 
size of the x-ray emitting region are still much 
larger than the light day timescales inferred for 
the TeV region. 

Early in 2008 Feb we detected an increase in 
the nuclear X-ray emission to a level a bit higher 
than it has ever been (figs. l8|9( ). This single Chan- 
dra observation was made during the TeV fiaring 
observed by MAGIC and VERITAS and the cor- 
responding values of fpy were f py~^ =+11. 8±0. 6 
and fpy- =-8.0±0.5. So far in 2008, HST-1 has 
been at a low level with only small changes in 
amplitude. Unfortunately we do not have a good 
estimate of X-ray timescales for the 2005 HESS 
event. Although we had a series of weekly obser- 
vations, the large fyp(nuclear) values could have 
been contaminated by rapid changes in HST-1, 
as indicated by fpy(HST-l) measurements derived 
from the readout streak photometry. If there is 
substantial TeV fiaring in 2009, we may be able 
to settle this question via an approved target of 
opportunity Chandra proposal which aims to find 
correlated X-ray /TeV behavior in the respective 
light curves. 

9. Summary 

We have found a quasi-periodic impulsive sig- 
nature in the brightening and dimming of HST-1 
in the X-rays. While this could be interpreted as 
a manifestation of past modulation of jet power, 
we suspect that it is rather a local oscillation of 
the process that converts bulk kinetic jet power to 
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Fig. 9. — X-ray lightcurves for the nucleus 
(squares), HST-1 (circles), knots D (triangles) and 
A (diamonds). 5% of the intensity of HST-1 has 
been subtracted from the nuclear values. We con- 
sider it highly probable that most / all of the knot D 
apparent variability is simply contamination from 
HST-1, and the slight shift to later time of the 
peak in 2005 is caused by the secondary response 
(release of trapped charge) of the HST-1 PSF. The 
secular decline of the knot A lightcurve is roughly 
consistent with the loss of effective area at low 
energies caused by contamination buildup on the 
ACIS filter. 

the internal energy of the emitting plasma. The 
fact that HST-1 lies on the northern edge of the 
cone defined by the VLBA jet (Paper IV) and that 
the cross section of HST-1 is less than 0.1% of the 
jet area (the cone of the VLBA jet has a diame- 
ter of ~ 90 mas at the distance of HST-1 whereas 
the effective diameter of the unresolved upstream 
end of HST-1 is ~ 2.5 mas), leads us to speculate 
that the time varying acceleration of electrons is 
related to a local instability. The ratio of <0.1% 
in areas is consistent with the ratio of power emit- 
ted by HST-1 to that believed to be the k inetic 
power of the jet ( Bicknell fc BegelmanI 1996h . 

The finding that the decay times of the lightcurves 
of HST-1 progressively lengthen moving from high 



to low frequencies suggests that simple expansion 
is not the primary energy loss mechanism for the 
relativistic electrons. If the X-ray decay time actu- 
ally reflects the synchrotron halflife of the highest 
energy electrons, we are not only able to estimate 
a value of the magnetic field which is independent 
of the usual equipartition assumptions, but also 
to provide an explanation of why the impulsive 
brightening is seen only at X-rays. That is be- 
cause the period of the oscillations is very similar 
to the X-ray decay time. The UV decay time is of 
order 10 times longer, so that an oscillating bright- 
ening would be smoothed out by the failure of the 
UV radiating electrons to lose their energy before 
the next brightening. The similarity of the decay 
time to the oscillatory brightening (both at X-ray 
frequencies) may be a coincidence (which makes 
the effect manifest), or it could be a causal com- 
ponent of the (as yet to be determined) instability 
mechanism. 

The X-ray variability timescale evidence sug- 
gests that the nucleus displays faster variability 
than does HST-1. This is circumstantial evidence 
for the hypothesis that the site of the flaring TeV 
emissions is the unresolved nucleus rather than 
HST-1. However, the shortest nuclear variability 
timescale we can measure from the Chandra data 
(< 20 days) is still significantly longer than the 
shortest TeV variability of M 87 reported by the 
HESS and MAGIC telescopes (1 - 2 days). 

The analogy we find useful in thinking about 
HST-1 is that of a river representing the underly- 
ing power flow of the jet. When something occurs 
to transfer some fraction of this power flow into a 
radiating plasma, we think of this as 'white water'. 
While knot A might be compared to a waterfall, 
HST-1 is more like a rock in the river. Thus the 
giant flare could have been caused by a change in 
local conditions (e.g. something moving into the 
river), with the resulting bits of white water car- 
ried downstream as the observed radio blobs. 

We acknowledge the Chandra Director's Office 
for approving our request to obtain a series of ob- 
servations with short intervals in 2007 Feb and 
Mar. Those data were of critical importance in 
distinguishing timescale differences between the 
nucleus and HST-1. A. S. Wilson (deceased) & 
W. Sparks were co-I's on the yearly Chandra pro- 
posals. We have benefited from discussions with 



11 



R. Blandford; with colleagues in the 'TeV com- 
munity', particularly H. Krawczynski, P. Collin, 
and M. Hui; with F. Primini on aspects of error 
propagation; and with B. Wargelin, G. Allen, and 
R. Edgar for advice on the subtle aspects of the 
ACIS CCD. We thank the anonymous referee for 
a useful critique. 

Work at SAO was supported by NASA grants 
G06-7112X, G07-8119X, and G08-9116X. C.C.C. 
acknowledges support from NRAO through a Jan- 
sky Postdoctoral Fellowship (2004-2007) and the 
NASA Postdoctoral Program at Goddard Space 
Flight Center, administered by Oak Ridge As- 
sociated Universities through a contract with 
NASA. L. S. acknowledges support by MEiN 
through the research project 1-P03D-003-29 in 
years 2005-2008. E.S.P. acknowledges support 
from the NASA LTSA program through grant 
NNX07AM17G. 

Facilities: VLA, HST, CXO (ACIS) 



12 



A. Pileup Problems 
A.l. Saturation 

Although the methods of § [5] recover most of the intensity of piled sources, there are on-board filters 
that keep events with energies > 17 keV or with certain bad grades from being telemetered to the ground. 
We suspect that these filters produced a significant decrease in our measurements when HST-1 was near its 
peak intensity, but we have not been able to verify this directly since the housekeeping files that provide the 
number of events 'dropped' because of amplitude and grade appear to be dominated by cosmic ray events. 

A.2. "Eat Thy Neighbor" 

When dealing with emission regions closer than an arcsec to each other, another pernicious effect of pileup 
comes from the detection algorithm. Whenever a candidate event is found, any other event within its 3x3 
pixel grid (for FAINT MODE) is considered to be a part of that event. Of the 2 (or more), the event with 
the most energy 'wins' and its position determines the reported location, and its energy is found from the 
sum of the charges within the 3x3 grid. HST-1 and the nucleus are separated by 0.86" or 1.75 pixels. Thus in 
the normal course of events, there will be occasional conflicts of this sort involving a nuclear count arriving 
within the same frame time as one from HST-1. As the intensity of HST-1 rises, this will happen more often, 
and many times the nucleus will 'win' because it has a harder spectrum than does HST-1. However, when 
pileup gets stronger, HST-1 will 'win' most of the time since almost all events in a frame time will consist of 
at least two photons, boosting the recorded charge so as to exceed the single photon from the nucleus. We 
are unaware of any quantitative estimates of this effect, but suspect it is not causing any serious problems 
for our analyses when HST-1 is < 4 keV/s. 

A. 3. Second Order EfTects of Pileup 

When there is negligible pileup, the mutual overlap of the core and HST-1 PSFs is « 5±2% for the 
rectangular regions used for photometry (fig. [T]). However, when pileup is significant, there are second order 
effects which seriously distort the PSF. Although we avoid some of these by summing energy instead of just 
counts, there are others which cannot be accommodated: e.g., the release of trapped charge during readout. 
Since the most obvious of these produces a secondary response displaced a few pixels from the PSF in the 
direction away from the readout buffer, the primary effect on adjacent jet features changes with the roll angle. 
Given the celestial position of M87, the roll angle is such that there are long periods of relatively constant 
roll angles at the beginning and end of the M87 viewing season (Nov. to Aug.) and a rather rapid change 
by close to 180° centered around a date late in March. Since the PA of the readout streak is similar to that 
of the jet at the beginning and end of the season, the primary second order effect produces contamination of 
jet features adjacent to HST-1 in the sense that prior to March the nucleus is badly affected whereas after 
the end of March it is knot D which suffers. This effect is evident in fig. [SI 

The fact that this sort of second order effect of pileup has not been calibrated (nor may be susceptible to 
calibration) means that changes in the measured intensity of the nucleus (and knot D) may not be intrinsic 
during the time that HST-1 was bright. For that reason, we have restricted usage of the core data to times 
when HST-1 was not too intense. 



B. X-ray Intensities for the Nucleus and HST-1 

In Paper III, we gave our measured intensities for HST-1 through 2005. Here we repeat these and also 
provide uncertainties and the values for the nucleus (already with 5% of HST-1 subtracted), both through 
the current season which ended in 2008 Aug. As described in Paper III, we can estimate the 'fudge factor', 
a, in the effective area by measuring the flux from fluxmaps when pileup is not serious. However, we plan 
to deal with spectral properties of the various features in a future paper so prefer to publish here only our 
directly determined intensities. 
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Table 1 

Chandra DATES AND X-RAY INTENSITIES OF THE NrCLEUS AND HST-1 





Observational Parameters 


Nucleus 


HST-1 


Epoch 


Date ObsID Livetinic 






label 


(sec) 


(kcV/s) (kcV/s) 


(kcV/s) (koV/s) 





2000 


Jul 30 


1808 


12845 


0, 


.300 


0, 


.005 


0. 


.247 


0, 


.005 


B 


2002 


Jan 16 


3085 


4889 


0, 


.644 


0, 


.014 


0. 


.734 


0, 


.013 


C 


2002 


Feb 12 


3084 


4655 


0, 


.575 


0, 


.013 


0. 


.588 


0, 


.012 


D 


2002 


Mar 30 


3086 


5089 


0, 


.536 


0, 


.013 


0. 


.576 


0, 


.012 


E 


2002 


Jun 08 


3087 


4973 


0, 


.443 


0, 


.012 


0. 


.798 


0, 


.014 


F 


2002 


Jul 24 


3088 


4708 


0, 


.494 


0, 


.013 


1, 


.096 


0, 


.017 


G 


2002 


Nov 17 


3975 


5287 


0, 


.649 


0, 


,014 


0. 


.799 


0, 


.013 


H 


2002 


Dec 29 


3976 


4792 


0, 


.652 


0. 


.014 


0. 


.653 


0. 


.013 


I 


2003 


Feb 04 


3977 


5276 


0, 


.598 


0, 


,013 


0. 


.645 


0, 


.012 


J 


2003 


Mar 09 


3978 


4852 


0, 


.777 


0, 


,016 


0. 


.872 


0, 


.015 


K 


2003 


Apr 14 


3979 


4492 


0, 


.583 


0, 


,014 


1. 


.071 


0, 


.017 


L 


2003 


May 18 


3980 


4788 


0, 


.378 


0, 


,011 


1. 


.022 


0, 


.016 


M 


2003 


Jul 03 


3981 


4677 


0, 


.372 


0. 


,011 


0. 


.859 


0. 


,015 


N 


2003 


Aug 08 


3982 


4841 


0, 


.304 


0. 


,010 


1. 


.214 


0. 


,018 





2003 


Nov 11 


4917 


5028 


0, 


.723 


0, 


,016 


2, 


.192 


0, 


,026 


P 


2003 


Dec 29 


4918 


4677 


0. 


.384 


0. 


,012 


2. 


.041 


0. 


,026 


Q 


2004 


Feb 12 


4919 


4703 


0, 


.515 


0. 


,016 


4. 


.079 


0. 


,04 




2004 


Mar 29 


4920 


5235 


















S 


2004 


May 13 


4921 


5251 


0, 


.396 


0. 


,013 


5. 


.358 


0. 


,05 


T 


2004 


Jun 23 


4922 


4543 


0, 


.228 


0. 


,013 


5. 


.323 


0. 


,05 


U 


2004 


Aug 05 


4923 


4633 


0, 


.264 


0, 


,014 


5. 


.636 


0, 


,055 


V 


2004 


Nov 26 


5737 


4237 


0, 


.740 


0. 


,022 


7. 


.494 


0. 


,07 


w 


2005 


Jan 24 


5738 


4666 


0, 


.416 


0. 


,023 


8. 


.316 


0. 


,08 


X 


2005 


Feb 14 


5739 


5154 


0, 


.439 


0, 


,024 


8. 


.785 


0, 


,08 


Ya 


2005 


Apr 22 


5740 


4699 


0, 


.254 


0. 


,019 


12. 


.417 


0. 


,11 


Yb 


2005 


Apr 28 


5744 


4699 


0, 


.505 


0. 


,022 


12. 


.167 


0. 


,11 


Yc 


2005 


May 04 


5745 


4705 


0, 


.435 


0, 


,021 


11, 


.798 


0, 


,108 


Yd 


2005 


May 13 


5746 


5142 


0, 


.557 


0. 


,022 


11. 


.555 


0. 


,10 


Ye 


2005 


May 22 


5747 


4701 


0, 


.597 


0. 


,023 


11. 


.243 


0. 


,10 


Yf 


2005 


May 30 


5748 


4699 


0, 


.268 


0, 


,018 


10. 


.665 


0, 


,10 


Yg 


2005 


Jun 03 


5741 


4698 


0, 


.357 


0. 


,019 


10. 


.432 


0. 


,095 


Yh 


2005 


Jun 21 


5742 


4703 


0, 


.252 


0. 


,018 


10. 


.270 


0. 


,094 


Yi 


2005 


Aug 06 


5743 


4672 


0, 


.165 


0, 


,013 


7. 


.098 


0, 


,067 


Yj 


2005 


Nov 29 


6299 


4655 


0, 


.514 


0. 


,016 


4. 


.032 


0. 


,042 


Yk 


2006 


Jan 04 


6300 


4660 


0, 


.612 


0. 


,017 


3. 


.833 


0. 


,040 


Yl 


2006 


Fob 19 


6301 


4337 


0, 


.861 


0, 


,020 


3. 


.535 


0, 


,040 


Ym 


2006 


Mar 30 


6302 


4701 


0, 


.407 


0. 


,014 


3. 


.786 


0. 


,040 


Yn 


2006 


May 21 


6303 


4699 


0, 


,377 


0. 


,013 


3. 


,161 


0. 


,035 


Yo 


2006 


Jun 19 


6304 


4677 


1, 


.117 


0. 


,021 


2. 


,557 


0. 


,031 
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Table 1 — Continued 



Observational Parameters Nucleus HST-1 



Epoch 




UDSIU 


Livetime 




<7 


i\ 


(7 


laoei 






(sec) 


[KS V / S J 


[Ke V / s j 


lyKe V / s j 


[KG V / S j 


Yp 


2006 Aug 02 


6305 


4653 


0.802 


0.018 


2.246 


0.028 


Yq 


2006 Nov 13 


7348 


4543 


0.682 


0.018 


4.042 


0.04 


Yr 


2007 Jan 04 


7349 


4685 


0.703 


0.018 


3.642 


0.039 


Ys 


2007 Feb 13 


7350 


4662 


0.823 


0.019 


3.516 


0.039 


Yt 


2007 Feb 15 


8510 


4701 


0.641 


0.017 


3.428 


0.038 


Yu 


2007 Feb 18 


8511 


4703 


0.561 


0.016 


3.515 


0.038 


Yv 


2007 Feb 21 


8512 


4703 


0.694 


0.017 


3.479 


0.038 


Yw 


2007 Feb 24 


8513 


4700 


0.717 


0.018 


3.472 


0.038 


Yx 


2007 Mar 12 


8514 


4471 


0.722 


0.018 


3.443 


0.039 


Yy 


2007 Mar 14 


8515 


4696 


0.651 


0.017 


3.533 


0.038 


Yz 


2007 Mar 19 


8516 


4679 


0.716 


0.018 


3.599 


0.039 


Za 


2007 Mar 22 


8517 


4674 


0.775 


0.019 


3.571 


0.039 


Zb 


2007 Mar 24 


7351 


4683 


0.728 


0.018 


3.396 


0.039 


Zc 


2007 May 15 


7352 


4588 


0.432 


0.014 


2.968 


0.034 


Zd 


2007 Jun 25 


7353 


4543 


0.372 


0.014 


3.445 


0.038 


Ze 


2007 Jul 31 


7354 


4707 


0.236 


0.011 


3.383 


0.037 


Zf 


2007 Nov 25 


8575 


4679 


0.610 


0.015 


1.376 


0.020 


Zg 


2008 Jan 05 


8576 


4692 


0.628 


0.015 


1.230 


0.019 


Zh 


2008 Feb 16 


8577 


4659 


1.480 


0.024 


1.305 


0.020 


Zi 


2008 Apr 01 


8578 


4706 


0.743 


0.017 


1.493 


0.021 


Zj 


2008 May 15 


8579 


4706 


0.565 


0.014 


1.398 


0.021 


Zk 


2008 Jun 24 


8580 


4705 


0.950 


0.019 


1.293 


0.020 


Zl 


2008 Aug 07 


8581 


4657 


0.431 


0.012 


0.947 


0.016 



'^The values of K are 'detector-based' observed intensities. They come from summing the 
energies of all events (from the 'evtl' file) within a rectangle of length 5 pixels transverse 
to the jet and 2.5 pixels along the jet (fig. [T]). No background subtraction was employed: 
when HST-1 is weak, the background (in counts) is of order 1%; for the highest intensity, the 
background is less than 0.3%. There is no correction for the build-up of ACIS contamination. 

'^Uncertainties are based on the raw counts: ^/N /N, hence la. 

Archival data from I Wilson fc Yantd (|2002l ). 



''This observation was taken in continuous clocking mode, so there is no 2D image available. 
All other observations had l/8th subarray ACIS-S3 chip only, and 0.4s frame time. 
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